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Figure 1: Sketch of the overall strategy to determine the A-parameter in Af = 3 QCD {left) as well as the 
renormalization group running to invariant quark mass M {right). 


1. Introduction 


In a mass-independent renormalization scheme for QCD, the renormalization group (RG) 
equations for the running coupling and quark mass read 



g^O 

-g^{bo + big^ + b2^ + ...) , 

(1.1) 

^^m{n) = T{g) 


—g^{do + dig^ + ••■)) 

( 1 . 2 ) 


with universal coefficients bo,bi,do and higher order ones, and dj^o, which are scheme- 
dependent. In such a scheme it is natural to first solve the RG equation for the coupling (1.1) 
and then for the mass (1.2), as the latter depends parametrically on the coupling only. Formally, 
the solutions can be exactly written down via renormalization group invariants (RGI), 


A = /I [bog^{lJ.)] e G(2 W(m)) exp | 

Mi = fTiiiiJ.) [ibog^il^)] exp |-^ 



■ 1 1 bi' 

p{g)^bog^ big 



i e {u,d,s,c,b,t}, 


(1.3) 

(1.4) 


which encode all information about the respective fundamental parameters of QCD. They are de¬ 
fined withouf relying on perturbation theory, and, while the Lambda-parameter is trivially scheme- 
dependent, the RGI mass is independent of the renormalization scheme. Accordingly, the RGI pa¬ 
rameters allow for easy conversion (at high but finite /i) between renormalized masses or couplings 
in different (massless) schemes. The only remnant and non-trivial dependence is on the number 
of quark flavours Nf. The RGIs {A,M,}a?, have been computed by the ALPHA collaboration in 
the past for Af = 0,2 QCD, cf. Refs. [1, 2, 3, 4, 5]. In line with the current Coordinated Lattice 
Simulation effort (CLS) [6], a new computation has been started to determine the fundamental pa¬ 
rameters {A,M,};Vf for QCD with Af = 3 dynamical flavours. The standard techniques involved 
here have been refined over the years. As intermediate mass-independent renormalization scheme 
the Schrddinger functional (SF) is used with non-perturbatively improved Wilson fermions. In the 
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mSF 

Lj a 

/3 

^SF 

^crit 

“fit 

1.1100 

6 

8.5403(55) 

0.1323361(12) 

1.1100(12) 

1.1100 

8 

8.7325(72) 

0.1321338(13) 

1.1100(15) 

1.1100 

12 

8.995(11) 

0.1318617(10) 

1.1100(24) 

1.4808 

6 

7.2618(28) 

0.1339337(13) 

1.4808(11) 

1.4808 

8 

7.4424(38) 

0.1336745(11) 

1.4808(15) 

1.4808 

12 

7.7299(89) 

0.13326299(69) 

1.4808^5) 

2.0120 

6 

6.2735(44) 

0.1355713(17) 

2.0120(32) 

2.0120 

8 

6.4680(51) 

0.1352362(15) 

2.0120(39) 

2.0120 

12 

6.7299(68) 

0.1347591(10) 

2.0120(49) 


Table 1: Simulation parameters and their known accuracy. 


SF the renormalization scale is given by the inverse box size /i = 1/L. Combined with a recursive 
finite-size scaling L/s in the continuum, typically with s = 2, one can reach a scale difference 
of O(IO^) after = 10 steps, and thus bridge the gap between typical large volume simulations 
(L ~ 6fm) — where the scale has been set — and the regime where perturbation theory can be 
safely applied. 


2. Revised strategy for Nf = 3 


In order to achieve a precision on as{Mz) that is better than the current PGD estimate, the 
ALPHA collaboration has elaborated a new strategy which minimizes statistical and systematic 
effects by combining the standard SF coupling scheme with the new gradient flow (GF) sheme in 
the Schrddinger functional [7, 8]. To this end one has to match the two schemes non-perturbatively 
at some intermediate scale /Xswi = 1 /Lswi, as indicated in the sketch shown in the left panel of Fig. 1 
and explained in more detail at last year’s Lattice conference [9]. A status report has been given 
this year [10]. In Fig. 1 the red crosses indicate values for the gradient flow coupling = gQp(L), 
L > Lswi, which is used to make contact between a hadronic scale in large volumes, made available 
by CLS, and the switching scale Ls^i- From there on the Schrddinger functional coupling = 
gsp(L), L < Lswi, is used (blue circles) to make contact with perturbation theory (dashed curve) at 
very high energy in order to extract the A-parameter in physical units according to eq. (1.1). 

As aforementioned, solving the RG equation for the strong coupling is a prerequesite for 
solving other RG equations such as the one for the mass which we want to consider now. The 
general strategy for computing RGI quark masses follows the decomposition [4, 5] 


M: = 


M 


/M(/rhad) 


X m,’(/rhad) , 


Mhad — 1/Thad , (2-1) 


where Lhad is some hadronic scale of 0(1 fm) and the total RG running factor for the mass in the 
continuum, M/m(jUhad), connects the renormalized current quark mass m;(fthad) to its RGI value. 
This factor does not depend on the individual quark flavour / e {u/d,s,c,b}, and in the Af = 2 
case, with A[M/hi] = 1.1%, it has been one of the dominant sources of errors in the determination 
of Ms [3] and Mb [11]. Accordingly, it is advisable to reduce this error systematically as much as 
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possible in the forthcoming determination. Due to the recursive step-scaling procedure the total 
RG factor decomposes into 


M 

iniUhad) 


M 


N 

X ]^ap(M/), 
1=1 


(^had/2') , ftpT = 2^l^had , (2-2) 


i.e., a product of step-scaling functions (SSF), ap{ui), and a part that can be safely evaluated in 
perturbation theory using eq. (1.4) after a sufficiently high scale /Tpt is reached. The SSF reads 


<yp{ui) = exp 



^(g) 

Hig) 




limrp(M,-,a/L) , 

a^O 


Zp{ui,a/L) 


Zp(go,2L/a) 

^p{go,L/a) 


i.e., it is obtained from its lattice SSF Zp along a line of constant physics (LCP) when a —)■ 0. This 
is achieved by implicitly prescribing a value u, to the renormalized coupling in a given scheme 
at vanishing quark mass niq. The relevant scale-dependent renormalization constant Zp is obtained 
from standard coiTelation functions (/p,/i) in the SF by the renormalization condition 

Zp{gQ,Lla) ^'^^^^ =C3(a/L), (2.4) 

V-^/l J m,=0 


employing the known tree-level normalization c^{a/L) with vanishing boundary gauge fields and 
(r /L,d,mq) = (1,0.5,0), cf. [4, 5]. For fhe Ihree-flavour compulation we are following as closely 
as possible fhe slrafegy for fhe running coupling, meaning fhaf we also swilch fhe scheme for selling 
up Ihe LCP above/below Lgwi- The compulation of Ihe lolal RG factor is easily extended to Ihis 
Iwo-scheme case (GF/SF coupling) and reads 


M _ M ^ m(AtpT) ^ m(/iswi) 
wi(Athad) Wi(AtPT) '^(Atswi) wi(Athad) ’ 

wilh 


AtPT = 2'^^pAiswi , Mswi = 2"^® Athad, (2.5a) 


m(AtpT) 


VsF 


nap(nP) 


Wi(Mhad) 


Vgf 




,?SF(2"swi/2 ) , 

(2.5b) 

fGF(2-had/2') . 

(2.5c) 


By doing so we direclly prolil from Ihe running coupling projecl [9, 10] as Ihe luning to vanishing 
quark mass and fixed coupling has already been done very accurately in terms of fhe bare parame¬ 
ters (j8, tq;rit) of Ihe Iheory. An example for Ihe parameters in use is given for Ihe SF coupling pari 
in Table 1. Fixing Ihe SF coupling to Ihe values in Ihe firsl column allows for individual Lja to 
use Ihe known paramelrizalion (j3) to determine Ihe relevanl value of j8 = 6/gQ, as given wilh 
error in column 3. Our knowledge of Ihe critical line where ntq = 0, subsequenlly gives 

Ihe critical hopping parameters as listed in column 4. Column 5 finally quotes Ihe value of (j8) 
al Ihe previously determined value of j8 and Ihus reflecls our knowledge of Ihe SF coupling. 


3. Preliminary results 

In Ihe Af = 3 running coupling projecl [9, 10], Ihe step-scaling function (SSF) for Ihe RG 
running of Ihe SF coupling in Ihe high-energy regime, asF(M^^) (2.012 > > 1.110), has been 
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Figure 2: Individual (prelimi¬ 

nary) continuum extrapolations of 
T,p{u^^,a/L) and its one-loop im¬ 
proved results 'e!'p\u^^, a/L) using all 
available lattices at fixed SF coupling 
{2.012,1.4808,1.110}. 
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Figure 3: Comparison of step-scaling functions ap(M^^) for Nf = 2 (red diamonds) with our preliminary 
results in the three-flavour theory (blue squares). The solid lines represent the perturbative running using the 
2-loop T-function and 3-loop j3-function in the 2- and 3-flavour theory, respectively. 


finished. Its determination requires field configurations with non-vanishing boundary conditions 
on the gauge field [12], while determining renormalization factors, such as Zp and thus Zp, usually 
proceeds with vanishing boundary gauge fields. To this end we have started production of N{ = 3 
field ensembles for L/a G (6,8,12} and their ILja counterparts. As for the step-scaling function 
asF this is done at 8 values of covering the aforementioned range. The bare parameters for 
these runs are summarized in Table 1 for the three values G {2.012,1.4808,1.110} at which 
we already have results on the computationally more expensive 2L/a = 24 lattices. We present 
preliminary results for the lattice step-scaling l.p{u,a/L) in Figure 2 (blue filled circles), together 
with its 1-loop improved version (red filled squares), defined by 


LP{u,a/L) 


Lp{u,a/L) 

1 -|- 5p{a/L) u 


(3.1) 


This effectively reduces the leading cutoff effects in the lattice step-scaling function 12p{u,a/L). 
The coefficients 8p{a/L) can be inferred from [13]. Although, it is not clear whether the rightmost 
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points for rp(M,l/6) are still in the region where (a/L)”>^ cutoff effects are negligible, cf. [5], 
for the time being we use all three data points and perform the continuum extrapolation using the 
ansatz rp(M,a/L) =A + B{a/L)^ to determine A = ap(u). We obtain 


u 

2.012 

1.4808 

1.110 

Gp{u) from Zp 

0.9173(11)(19) 

0.9426(8)(17) 

0.9584(9)(17) 

Gp{u) from Zp ^ 

0.9160(11)(8) 

0.9416(8)(9) 

0.9576(9)(9) 


where the first error is statistical and the is the difference to the result of a 2-pt weighted average 
which neglects the L/a = 6 data. Our results at fixed u agree well at the one-sigma level. As seen 
in Figure 2, cutoff effects are reduced in the continuum extrapolation for Zp as to be expected. 
Accordingly, we can take the values in the third row as our preliminary results for ap(M^^) and 
note that at present the statistical, 5stat, and systematic eiTor, 5sys, are < l%o. To get an impression 
about the quality of these results, we compare them to results obtained for the mass SSF in the 
two-flavour case. To this end we reproduce Fig.3 of [5] in the left panel of Figure 3 above. We 
add to the Af = 2 data (red diamonds) our present results (blue squares) and note that the inner 
error bar is the purely statistical eiTor, and the outer error bars include the linearly added systematic 
errors. The latter have been determined as discussed earlier. The two solid lines represent the 
perturbatively known SSF ap^(M) to the highest available order in PT, cf. [5], for the Af = 2,3 
flavour cases respectively. The Af = 3 data falls in line with the PT SSF at small couplings in the 
SF scheme, best seen in the right panel of the same figure. We furthermore indicate the distribution 
of the couplings up = g|p(Lniax/2') as they have been in the Af = 2 case [3]. 

The fact that the new data scatters less about the PT behaviour at small couplings compared 
to Af = 2, is most likely a result of the intensified and more systematic effort to set up the LCP [9] 
which has become possible due to algorithmic and computational advancements during the last 
10 years. Especially at the switching scale = 2.012), also the total error could be reduced 
significantly. With 5 additional, uncorrelated continuum data points to be added to this picture 
soon, uf^ G (1.110,2.012), we will be in the comfortable position to fit the non-perturbative SSF 
to much higher accuracy than in the past. However, to give an upper bound on the estimate of 
the total error of the full Af = 3 RG running mass factor right now, we can naively assume, using 
A[ap] = 5stat + 5sys < 2%o as discussed above, that 


A 


M 

/M(/rhad) 


N 

]^A[ap(«,•)] ~ Vn-27oo ~ 


0.49% 

0.63% 


for A = 6 (as for Af = 2) 
for A = 10 


(3.2) 


This means that the total error compared to Af = 2 will be reduced by at least a factor of 2. Note 
that we have A = Asp + Agf according to (2.5). The main assumption entering this upper bound 
is that we control the error on ap(M,) as good in the GF scheme running part at lower energies as 
we have reported here for the SF scheme running mass part. Needless to say that the naive bound 
in (3.2) can be further reduced by increasing the statistics on the lattice SSF to reduce the statictical 
error even further. But as the systematic error is of the same size as the statistical, it would be 
better to add another lattice spacing such as L/a = 10 or 16 in order to check and further reduce 
the systematic error in the continuum extrapolation. 
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4. Summary 


We have presented first results towards a computation of the full RG factor for the running 
mass in three-flavour QCD by the ALPHA collaboration. This is a necessary ingredient to deter¬ 
mine renormalization group invariant quark masses from Af = 3 large-volume ensembles provided 
by CLS. Systematic errors due to large scale differences are controlled by following a traditional 
finite-size scaling approach, using the Schrddinger functional as intermediate finite-volume renor¬ 
malization scheme to non-perturbatively connect low- and high-energy regimes. 

A gain in the overall precision is achieved by (a) exploiting very precise results from the 
ongoing Af = 3 running coupling project, incorporating a non-perturbative scheme switch at in¬ 
termediate energies, and (b) via computational advances in algorithms and hardware over the last 
decade. Both together will allow us to reach an unprecedented precision in the non-perturbative 
step-scaling function of the mass, provided that equally good results are obtained for the running at 
lower energies. However, to finally quote an RGI quark mass in physical units, one needs to incor¬ 
porate the scale-setting procedure which introduces yet another uncertainty and modifies (2.1) to 


Mi = 


M 

— -- X hm 

Wl(Athad) 


amiil^had) 

^/had 


^ /had ) 


Mhad — 1 /f^had • (4.1) 


Here /had is any physical quantity that has been used to set the scale on the CLS ensembles. 

Together with improvements in the scale-setting procedure w.r.t. the two-flavour case, we can 
remain optimistic to reach for instance an overall precision less than 1 % on Mg itself. 
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